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Abstract 

In this paper plane elastic curves are revisited from a viewpoint 
that emphasizes curvature properties of these curves. The family of 
elastic curves is considered in dependence of a tension parameter a and 
the squared global curvature maximum It is shown that for any 
elastic curve is bigger than the tension parameter a. A curvature 
analysis of the fundamental forms of the elastic curves is presented. 
A formula is established that gives the maximum turning angle of an 
elastica as a function depending on and a. Finally, it is shown 
that an elastic curve can be represented as a linear combination of 
its curvature, arc length and energy function and that any curve with 
this property is an elastica. 



1 Introduction 



The search for Smoothing Algorithms in CAGD has induced research on 
curves and surfaces which minimize functionals with geometrical or physical 
meaning. The classical example are elastica which describe the shape of 
elastic materials. Plane elastica can be defined as the extremals (critical 
points) of the variational problem 




+ <7 ds — ♦ 



min 



where k denotes the curvature, s the arc length of a plane curve x and a is 
constant. The total length L of x is considered to be variable. The set of 
comparison curves of the problem is the set of all C°° curves in the plane 
with fixed endpoints and fixed tangent directions at these points. 
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Lee and Forsythe found (see [6]) that the curvature function of an elastica 
satisfies the differential equation 

k"(s) + \k 3 (s) + Kk(s) = 0, KeR. (1) 

We begin this paper with a short derivation of this result and a proof that this 
equation is also sufficient for a curve to be an elastica. We will then show that 
for any elastica the square of the curvature function has a global maximum 
which is bigger than the tension parameter a. The fundamental forms of 
elastic curves first described by Euler are characterized by the values of the 
parameters k ^ and o. A curvature analysis of these curves is presented. 

In section 3 various formulas expressing properties of elastica in terms of 
and a are given. This includes a formula that gives the maximum turning 
angle of an elastica as a function depending on k ^ and a. The explicit 
formula for the maximum turning angle allows to determine an upper bound 
for in the case that an inflectional elastica has a turning angle bigger 
than n. As the main result of this section we show that a plane elastica has a 
representation as a linear combination of its curvature, arc length and energy 
function E(s) = / 0 J k 2 (s) ds. This property characterizes the class of plane 
elastica and can be used for speeding up the computation of interpolating 
elastica. 



2 The curvature of plane elastica 

Given two points P, Q € R 2 and two unit vectors V € T p R 2 , W € TqR 7 
then M denotes the set 

M := {x : [0, L] — ► R 2 : L € R + , x € C°°[ 0, L), |i'(s)| = 1 for s € [0, L], 

x(0) = P, x(L) = Q , x'(0) = V, x'(L) = W). 

We consider the problem of minimizing the functional 

rL 

E + crL := / k 2 (s) + a ds 

Jo 

among all curves of M where k denotes the curvature and s the arc length 
parameter of a plane curve x. The integral E = Jq k 2 (s) ds is called the en- 
ergy of x on [0, L). a denotes a constant which we call the tension parameter. 
The total length L is considered to be variable. As the functional considered 
is translation and rotation invariant we assume in this section that P = (0,0) 
and V = (1,0). 
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For x € M the tangent vector T of x is given by 
T(s ) = ( cos('P(s)), sin(ty(s)) ) 



where the function with 'P(s) := / 0 * k(s ) ds gives the turning angle of x. 
Using this notation the variational problem can be written as 



f L 

min / ^(s) + a ds 

*eC°°(0,L]./o 



under the constraint 



fra)*-* 

where the admissible functions 6 C°°[0, L\ are subject to the constraints: 

ty(0) = 0, 'I'(L) = ip with ( cosip , sin^’) = IV. 

According to Lagrange’s multiplier rule the differential equation of this 
problem is the Euler equation for the integrand 

F(^, V'js) = tf' 2 (s) + a + Acos'P(s) + /zsin 'I'(s), 

for some constants A,/z £ R, i.e. 



tf"(s) = — — sin 'l'(s) + — • cos 'I'(s). 

(see e.g. in [3], [6]). By introducing the constants 

A = 2 a cos <p /z = 2a sin (p 



( 2 ) 



we rewrite the Euler equation in the form 

ty" = n’ = —a sin(^ — <P). (3) 

Multiplying (3) by 2'P / and integrating yields 

k 2 = 2a cos('P — <P) + A (4) 

where A € R denotes an integration constant. Equation (4) has been used 
to define elastic curves in the plane in the classic literature (see [7]). 

In order to determine the constant A in terms of the tension parameter 
< 7 , we consider the boundary condition 

F(nL), *'(Z,),Z.) - ^;F(*(i), *'( 5 ), %¥'(/,) = 0 
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that must be satisfied by the extremal. This condition is implied by the fact 
that the total length L of the curve is variable in the variation (see e.g. [2], 
p. 571). Thus, 

^(L) = k 2 (L) = A cos 'I'(L) + p sin V(L) + <r. (5) 

Comparing (5) with (4) shows that 

A — a. 

Therefore we give the following definition of an elastic curve under tension. 

Definition 1 An arc length parametrized plane curve x with curvature func- 
tion k is called elastica (or elastic curve) with tension parameter a, if for 
some a, tp 6 R 

K 2 = 2a cos('P — <f>) + a (6) 

where denotes the function ty(s) := fo K,(s)ds. 

Lee and Forsythe showed in [6] that the Euler equation (2) implies the 
differential equation (1) for the curvature function k. We will now show that 
(1) is in fact equivalent to (6). 

Theorem 2 If k € C 2 (R) and 'I'(s) := f£ n(s)ds, then 

k 2 = 2a cos('P — 4>) + a 

holds for some constants a, ip and a if and only if 

k" + - ^a/c = 0. (7) 

Proof: (i) Differentiating (6) yields the Euler equation in the form of 
(3). Differentiating the Euler equation gives 

k" — —a cos('I' — 4>)k. 

Substituting the term cos(^ — y>) in the above equation according to (6) we 
obtain (7). 

(ii) Differential equation (7) can be integrated to the first oder differential 
equation 



(k') 2 = C — (l/4)a 2 — (l/4)(« 2 — a) 2 . 
Note, that for a real solution of (8) it is necessary that 

a 2 := C — (l/4)a 2 > 0 



( 8 ) 
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and 



(" a ~ *) 
4a 2 



< 1 . 



The function 6 defined by 



6:= (-irarccos(^— -) 

2\a\ 

with n = 1 for «' > 0 and n = 2 for k' < 0 obviously obeys the relation 

k 2 = 2\a \ cos 6 + a. (9) 

It remains to show that 6 is an integral of n. 

(8) together with the definition of a 2 and (9) implies 

(k') 2 = a 2 sin 2 6. 

Since sin 6 is positive (negative) if k' is negative (positive), we obtain 

k = — |o| sin 6. (10) 

Differentiating (9) and substituting k' according to (10) yields 

& = K. 



Since and 6 are both integrals of k there is a constant <f> such that 
$ = 6 + 4>. For 

k 2 = 2|a| cos(^ — <f>) + a 

holds because of (9). 

□ 



The squared curvature function k 2 of an elastic curve has a global maxi- 
mum, even if « 2 is extended to the whole real line. 

Lemma 3 If k is a solution of (7) on R, then k 2 has a global maximum. 

Proof: Since the function (1/2)(k 3 — crn) is of the class C J (R) the solu- 
tions of (7) can be extended to the whole real line. 

To show that k 2 has a global maximum we observe first that (3) implies 
the existence of a local extremum of «. This is because the assumption 
k'(s) = *"{s) ^ 0 for all s 6 R means that is convex or concave and 
therefore unbounded while according to (6) k' has no zeros only if $ is 
bounded. Furthermore from (3) and (4) it is obvious that any local extremum 
of k is a global extremum of k 2 . 
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We assume now that k 2 has no global maximum on R. In this situation 
k' has exactly one zero s m j n on R and k 2 takes its global minimum in this 
point. Note that k is non-zero for any s ^ s mtn because a zero of k at a 
point s ^ s m)n would imply that K(s mtn ) = 0 and k'(s) = 0 for some point s 
between s and s m)n . Therefore ( k 2 )' has no zero besides s m)n and k 2 is strictly 
monotone increasing on the right of s mtn . Hence n is monotone increasing 
resp. decreasing on the right of s min if k has positive resp. negative values 
on the right of s mtn . The formula 

y(s) = J K(s)ds + y{s min ) 

yields that 'I' is in any case unbounded. This is a contradiction to the as- 
sumption because (4) implies the existence of global maxima for « 2 if 'J' is 
unbounded. 

□ 



We now give the relation between the global maximum k ^ of re 2 and the 
tension parameter c and express k in terms of elliptic functions. 



Theorem 4 Let k € C 2 (R) be a solution of the differential equation (7) with 
a global maximum / 0 o/ k 2 on R. Then the following statements hold: 

(0 K m> a 

(ii) k has a zero if and only if >2 a . In this case k is given by 



k(s) ~ K m cn(^/(^ - c)/2(s - s m ) | k 2 ) 



( 11 ) 



with the parameter 



k 2 = 



2(4> - a) 

(iii) 2 c — n 2 m < k 2 < for n 2 m < 2c. In this case k is given by 
k(s) = K m dn(*c m (s - s m )/2 | p). 



( 12 ) 



Proof: As a global maximum point s m of « 2 is also a zero of k' (8) implies 



C = — k! — —ck 2 . 
4 m 2 



Hence (8) takes the form 



(k ') 2 = 7 («m - K 2 )(k 2 + «rn - 2<7 )- 



(13) 
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As all quantities in (13) are real (i) follows from the fact that the term 
( K m ” * 2 ) ’ s always non-negative and therefore the term (/c 2 + — 2 <r) has 

also to be non-negative for any k 2 . 

Using the same argument we observe that k 2 (s) = 0 for some s implies 
that > 2 a while for K 2 m < 2<r « 2 has to be greater or equal 2<x — k^. 

To express « in terms of elliptic functions we procede as follows. In the 
case that > 2<r we substitute z 2 = — k 2 )/«^ in (13) and obtain for z 

the differential equation 



<00 (i4) 

where 




(14) is the differential equation of Jacobi’s function sn=sn(u) for the argu- 
ment u = — <r)s (see [l],p. 114). We therefore obtain 



z 2 (s) = sn 2 ( v /(k^ 1 - <r)/2(s - s m ) \ k 2 ). 

The relation k 2 = k^(1 — z 2 ) then gives 

« 2 ( 5 ) = “ <t )/ 2 ( s ~ s m) I * 2 ) (15) 

(see [8],p.l6). Since k is differentiable, (15) implies (11). 

In the case that K 2 m < 2<r we substitute z 2 = ( k 2 m — k 2 )/2(k£, — a) in (14) 
and obtain for z the differential equation 

(*T = j«l(i - **)(i - iV) 



where 

Hence 

The relation 
yields 



2 1 2 (K 2 m - a) 



1 ~k 2 = 



z 2 (s) = sn 2 (« m (s - s m )/2 1 1 2 ). 
K 2 = K 2 m (\-l 2 Z 2 ) 



« 2 ( 5 ) = «m dn2 («m(5 ~ S m )/2 | l 2 ) 

(see [8],p.l6) which implies (13) again because of the differentiability of «. 
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Theorem 2 implies that the curvature function of an elastic curve extends 
to a periodic function on R. In the case that k ^ > 2 o this periodic function is 
symmetric with respect to any zero of its derivative and antisymmetric with 
respect to any zero. According to Love [7] this situation is called inflectional 
because the extension of x has turning points. In the case that k ^ < 2cr the 
extension of k has no zeros but is still symmetric with respect to any zero 
of its derivative. The local extrema in this case are K m and ±(2cr — /c^). 
The situation is illustrated in figure 1 - figure 8 where for « m = 1 various 
curvature functions and the corresponding elastica are shown. If only non- 
negative tension values are considered a is an element of [0,1]. With a 
increasing from 0 to 1 the curvature function changes continuously from a 
lemniscate function {a = 0) to a constant (<r = 1). A classification of the 
different forms of elastica has been given by Euler (see [5], [7]). 
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Figure 1. a = 0; k(s) = coElemn(/c m (a — s m )/2). 





Figure 2. a = 0.3; k ( s ) = ^ cn( x /(/c2 i - a)/2(s - s m ) | A: 2 ). 




Figure 4. c — 0.394757217; #c(s) as in figure 3. 



9 



10 * 





Figure 3. a- 0.49; «(s) = /c m cn^/c^ - cr)/2(s - £ m ) | it 2 ). 
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Figure A. a - 0.5; k(s) = Bech(/c tn (5 - s m )/2). 





Figure 5 m cr — 0 . 62 j /c(.s) — Kir. dn(/t fn (s 5 m )/2 | 



4l 

3 

2 ■ 

,1 

0 ■ 

.1 - 
.21 
-3 * 
u - 

**13 *10 -3 0 3 10 13 



3 
2 
1 
0 
*1 
-2 

.3 .1 0 1 3 s 4 




Figure 6. a = 1; k ( s ) = 
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3 A characterizing representation theorem 

First we summarize the main formulas for elastic curves. In order to avoid 
repetitions we use the following convention: if a curve x : [0,L] — > R 2 is an 
elastica , then k denotes the (analytic) curvature function of x, k denotes the 
extension of k on R, the global maximum of k 2 and s m is a number with 
K 2 (s m ) = Furthermore : [0, L] — ► R is the function 'I'(s) := / 0 * «(s)ds, 

* denotes the extension of 'I' on R and <p := ^(s,,,). The energy E of x is 
the function E(s ) = / 0 * k 2 (s) ds. For convenience we also use the notation 
C(tp) = (cos sin 0). i? is the angle in [0, 27r[ with x'(0) = C(t?). 

Theorem 5 For a plane elastic curve x with tension parameter a the fol- 
lowing relations hold: 

(i) k'(s) = -±(K 2 m - a) sin(^ - p), 

(ii) /c 2 (s) = ( k 2 m - a) cos(^(s) - p) + a, 

(iii) /c(s) - k(0) = “(*4 -a)< C'(<p + d),x{s) - x(0) >. 

(iv) E(s) = ( n 2 m -a) < C(<p + i?),x(s) - x(0) > +(ts, 

Proof: Without loss of generality we can assume that a in (4) is positive. 
From lemma 3 we know that k 2 has a global maximum k 2 m . Since for positive 
a the maximum occurs iff cos('I , (s m ) — <f>) = 1, it follows that 

2a = - o 



4> = $( 5 m ) 4 - 2kn . 

Therefore (3) and (4) imply (i) and (ii). To obtain (iii) we rewrite (i) as 

k '( s ) = “(«£, — <7)(sin( , I' + i?) cos(v? + i?) — cos('P + ) sin(v? + t?)) 

and integrate using 

x'(s) = (cos('P(s)),sin('I'(s))). 

To verify (iv) we write (ii) in the form 

/c 2 (s) = (k£j — <r)(cos(^ + t?) cos(p + i?) + sin(^ + i?) sin(<^> + t?)) + a 
and integrate. 

□ 
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Since the elliptic function cn in Theorem 3 has the symmetry properties 
of a sine wave while dn is positive, the turning angle of an elastic curve is 
bounded in the inflectional case a < and unbounded otherwise. Formula 
(ii) of Theorem 3 can be used to determine the maximum turning angle of 
an inflectional elastic curve. 



Corollary 6 The maximum turning angle ^ max := max <€ [ <)<+ r] | ^ (s) | of an 
elastica x with period T and tension parameter a < is given by 



Vmax = 2 arccos( 



O - nl 






(16) 



If a > ^ 0 then 

for any t € R. 



\V{t + T) - $(<)| = 2* 



Proof: Let s 0 be a zero of k. Then it follows from the symmetry prop- 
erties of k in the inflectional ca se that 

rto+T/4 

= 21/ k(s) ds\ = 2|'I'(s 0 + T/4) - ¥(a)|. 

J 30 

With (ii) we get 

AC 2 — (J —(7 

V max = 2| arccos(~Y ) + v? - arccos( — ) - <p\ 

K m- a K m ~ ° 

and therefore (16). 

In the non-inflectional case k is a positive periodic with the period T = 
4 A'/ 1 /c m | where K denotes the complete elliptic integral of the first kind. For 
t = s m one obtains 

|*(s m + T) — 'P (s m ) | = 4 / dn(u|/ 2 ) du = 4 arcsin(sn(A')) = 2n. 

Jo 

□ 



For an inflectional elastica x formula (16) implies an upper bound for n 2 m 
if the absolute value of the oriented angle between x'(0) and x'(L) is bigger 
than 7r. 

Corollary 7 For an inflectional elastica x with tension parameter o and 
ip := |^(L)| > tt : 

2 < — < ] !_. 

a cos(V , /2) 
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Proof: It follows from (16) that a turning angle bigger than 7r can only 
happen for positive a. The inflectional nature of x implies therefore the left 
inequality. The right inequality follows from the fact that ip has to be less 
or equal 'Pmax which is given by (16). 

□ 



As the main result of this paragraph we show that an elastica can be 
represented as a linear combination of its curvature, energy and arc length. 

Theorem 8 If x is an elastica with tension parameter a ^ then 

x(s \ _ 1 ( sin(p + tf) cos(^ + t?)\ /2 (k(s)-k(0))\ , QW 

K m - a ' - cos(y? + 7?) sin (y? + t?)/\ E(s) - as ) 

where cosy? = (/Cq — «t)/(k^ — a) and siny? = 2 Ko/(«in — a )- 

Proof: A plane curve x € C 00 ^, L ] has a representation 

x(s) — x(0) = / C(^(s) — y? + y? + i?) ds 

Jo 

= C(y? + i?) / cos('I / (s) — y?) ds 

Jo 

-f C"(y? + i?)/ sin('I , (s) — y?) ds . 

Jo 

Applying (iii) and (iv) of Theorem 3 yields 



x(s) - x(0) = -j- (C(v» + 1 9)(E(s) - as) - C'(<p + t?)2(«(s) - k(0))) 

K m ~ ° 

which is equivalent to (17). 

The formulas for follow from (i) and (ii) for s = 0. 



□ 



Formula (17) provides an explicit representation of an elastica in terms 
of its curvature function if n ^ ^ a. Note that in the case K 2 m = a the elastic 
curve is a circle of radius l/|« m | as shown in figure 8. 

Finally we show that elastica are essentially the only curves in the plane 
that have a representation of the form (17). 

Theorem 9 Let k be an arbitrary C 2 (R) function with a global maximum 
K 2 m of K 2 , a a real number smaller than n 2 m and E(s) := / 0 3 n 2 (s)ds. An arc 
length parametrized curve x where x(s) is given by (17) is an elastica with 
curvature function n and tension parameter a . 
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Proof: We assume the case + t? = 7r/2 which can always be achieved 
by applying a rotation to x. Then 

x' = — (2k', k 2 — cr), x" = —r~ (2k", 2 kk). 

* 1 -°- 

x is arc length parametrized if and only if x' = 1, i.e. 

(2 «') 2 = (4 - CT) 2 - (« 2 - <t) 2 . (18) 

Differentiating (18) yields 



2 k k" + k k(k 2 — <t) = 0. 

Therefore either k satisfies (7) or k is constant. For an arc length 
parametrized curve the determinant [x', x"] is curvature. Here we have 

[*',*"] = ^ - «"(« ! - <*))• (19) 

If k is constant the curvature of x is zero, hence x is a trivial elastic curve. 
If k satisfies (7) we substitute in (19) k" according to (7) and k' 2 according 
to (18). This yields 

r / 

[x , x J = AC. 

□ 

The representation formula (17) of an elastic curve is extremely useful for 
the computation of interpolating elastica. One reason for this is that (17) 
involves no trigonometric functions and fewer integrations than the standard 
representation based on the formula x' = (cos ty, sin ty). Therefore (17) is 
less expensive to evaluate. Furthermore (17) can be used to find piecewise 
polynomial approximations of elastic curves based on a spline approximation 
of the curvature function. The author has established polynomial splines 
which approximate the curvature functions of plane elastica with high ac- 
curacy. Using these piecew'ise polynomial curvature functions together with 
the new representation (17) one obtains polynomial spline approximations of 
the elastica itself. These approximations will be discussed in [4]. 
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